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Integrins mediate cell adhesion and are essential receptors for the develop- 
ment and functioning of multicellular organisms. Integrin activation is known 
to require both ligand and talin binding and to correlate with cluster forma- 
tion but the activation mechanism and precise roles of these processes are not 
yet resolved. Here mathematical modeling, with known experimental param- 
eters, is used to show that the binding of a stabilizing factor, such as talin, 
is alone insufficient to enable ligand-dependent integrin activation for all ob- 
served conditions; an additional positive feedback is required. 
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Introduction 



Integrins, large membrane- spanning heterodime ric proteins, were so named for their ability to 
link the extracellular and intracellular skeletons ITamkun et al., 1986*| . As an important class of 
cell adhesion receptors they participate in a wide-range of biological i nteractions, including de - 
velopment, tissue repair, angiogenesis, inflammation and hemostasis HHorwitz & Webb, 2003 1 . 
Cell adhesion and detachment as well as controlled actin polymerisation inside the cell are of 
particular importance in cell migration. The speed of cell movement depends on the density 
of integrins and ligands as well as their affinity of binding BPalecek et al, \99l\ . Integrins are 
key components of focal adhesions, dynamic multi-protein c omplexes that are invo lved in the 
regulation of cell adhesion and migration. Focal adhesions [Za mir & Geiger, 200 1| provide a 
physical link between integrins and the actin cytoskeleton as well as sites for signal transduction 
into the cell interior. Information about identified interactions and players in these complexes 
is ever-increasing but our overall understanding of how the ensemble works remains relatively 
poor. 

According to the current model of integrin activation, lig and binding sh ifts the equilibrium 
between different integrin conformations to the active one WHynes, 2002 1. The two extreme 
conformations of this allosteric protein are a bent or 'closed' conformation which represents 
the low affinity state for ligand and an 'open' conformation that will bind with high affinity 
to ligand. Conformational changes in the extracellular domain affect the cytoplasmic tails, 
which are separated in the open conformation but not in the closed. Separation of the cytoplas- 
mic domains promotes their interaction with cytoskeletal and signal transduction molecules, 
and thus the activation of integrins and downstream signaling. The conformational equilib- 
rium can be influenced both by ligand binding to the extracellular domain (outside-in signal- 
ing) and by binding of cytoplasmic proteins to the separated cytoplasmic domains (inside-out 
signaling). As well as changes in affinity induced by structural changes, integrins can also 
modulate their avidity by cl ustering, thus changing the valency of their interactions with ligand 
[ Carman & Springer, 2003| . 

Previous theore tical studies have addressed the mechanism of integrin clustering [Ward & Hammer, 19941 



Irvine et al. , 2002] but they have not yet explored whether ligand binding is sufficient or if other 



factors such as clustering are necessary for integrin activation. In principle, ligand-dependent 
integrin activation can be enabled by any of the four different mechanisms depicted in Figure 
El A) Extracellular ligand alone is sufficient to activate integrins and other factors are only im- 
portant for downstream signaling (Fig. EK); B) Ligand-dependent integrin activation requires 
the binding of an intracellular stabilizing factor that interacts with and stabilizes the on state 
(Fig. |2j3). C) A positive feed back, provid ed either by integrin-ligand pairs themselves or by 
self-interactions, is important [ |Li et al., 20031 ) (Fig. |2f!); D) A positive feedback that increases 
the activity of an intracellular stabilizing factor is required (Fig. EJ3). The latter would need to 
be based on a larger intracellular signaling network. 

Mathematical modeling is used in the following to evaluate the physiological potential of 
these mechanisms. The analysis suggests that, at least for some physiological conditions, bind- 
ing of ligand (even if it is preclustered) and/or a stabilizing factor (Fig. EK,B) are insufficient for 
ligand-dependent integrin activation; a positive feedback (Fig. EJl!,D) is required which results 
in the further stabilization of the active conformation. Integrin clustering is a likely consequence 
of such a positive feedback. 
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A model for integrin activation 



Overly simplified, integrins can be taken to exist in one of two states, a closed low affinity or an 
open high affinity conformation. In the absence of ligand the equilibrium is biased towards the 
inactive, low affinity conformation UTadokoro et ah, 2003| . Binding of ligand (L) stabilizes the 
active, high affinity conformation. Integrin activation can thus be measured as the fraction of in- 
tegrins bound by ligand. The analysis can be greatly simplified if the fraction of integrins bound 
to ligand is determined indirectly from the fraction of ligand bound by integrins. In this case, 
integrins in the closed, low affinity conformation do not have to be considered explicitly. We 
will assume that there is a constant number of unbound integrins in the high affinity conforma- 
tion since any integrin that binds to ligand can be replaced rapidly from the integrins in the low 
affinity conformation by conformational changes. While ligand binding will eventually deplete 
the pool of unbound integrins, this does not need to be included in an analysis of mechanisms 
for integrin signal initiation: if integrin-ligand binding occurs to a level that integrin depletion 
becomes relevant the activation mechanism can be considered successful independently of the 
exact final number of integrin-ligand complexes. 

Under physiological conditions, the ligand is expected to be sufficiently dense that spatial 
details and diffusion constraints can be neglected in the initial steps of focal adhesion assembly. 
Transport of integrins from parts of the cell that are not in contact with the substrate is not 
relevant for the initial activation of cell adhesion signaling cascades and is not considered in 
this model. 

The model follows the binding of single integrins to ligands; the probability of integrin- 
ligand complex formation can be calculated by solving the appropriate Master equation. The 
Master equation is a gain-loss equation for the probabilities of the separate states with n integrin- 
ligand complexes [van Kampen, 1992| , from which the macroscopic kinetic equation can be 
derived when the density fluctuations are negligible. In order to decide which of the four mech- 
anisms depicted in Figure|2]is sufficient for ligand-dependent integrin activation, the fraction of 
ligand bound by integrin is determined for each case. 



Ligand engagement is not sufficient for integrin activation 

According to the first mechanism (Fig. EKX ligand-dependent integrin activation does not re- 
quire any further supporting processes and interactions. The Master equation for the formation 
of i integrin-ligand complexes is thus given as 

Pi = k-(i + l)Pi+i + k+(L - i + - (k-i + k + (L - i))pi (1) 

with pi = if i £ [0, L], where pi is the probability that % (open) integrins are bound to a 
ligand with L sites. k + , k- refer to the on- and off-rate of the integrins. The on-rates include the 
density of integrins which are taken to be constant (see previous section). The first two terms in 
dU represent the gains of state i due to the transition from other states, that is due to unbinding 
of an integrin from a complex with i + 1 integrins (first term) or the binding to a complex with 
i — 1 integrins (second term). The bracketed term is the loss due to transitions from i into other 
states either by unbinding or binding of an integrin. 

For the steady state (j>i = 0) we obtain from ^ 




4 



The expectation value for the number of integrins bound to a given ligand with L binding 
sites is given by < >= J2i=o iPi- Integrin activation therefore requires ^> p (i > 0) and 
thus -jr- 3> 1 (Fig. 13), which corresponds to a high on- and low off-rate. This is because the 

probability of i bound integrins is proportional to the ith power of ^- and only for > 1 this 
probability increases with increasing i, that is with increasing numbers of integrins bound. 

k_ ~ 1 s _1 has been established experimentally [Vitte et al., 2004| ; the 2-dimensional on- 
rate is more difficult to determine and has to be calculated from the 2-dimensional dissociation 
constant which itself is difficult to measure and often obtained via its 3-dimensional counterpart 
(Kff ~6x 1(T 8 M for the high affinity conformation paull etal, 1993||Suehiro etal, 1997| 



and K"d > 1/iM for the low affinity conformation UFaull et al . , 1 993 ]). It is generally assumed 



that the two dissociation constants are linearly related s uch that K 2 ^ = rjK ^. The value of the 
conversion factor rj is, however, still a matter of debate [ |Moy et a l. , 1 999 1 , since theoretical es- 
timates are not in complete agreement with latest experimental results. According to theoretical 



estimates for a protein complex that spans about 20 nm r\ = 1.2 x 10 19 ^tt [Bell, 1978||. Studies 



employing protein-coated beads report r/ = 8 x 10 24 ^|p [Moy et al., 1999| and rj = 10 



122 



m 2 M 

[ |Kuo & Lauffenburger, 1993 j. The lower r\ value is corroborated by measurements of the 2- 
dimensional dissociati on constant for the L FA-3/CD2 adhesion pair that were carried out using 
cells instead of beads HDustin et al., 1996) . While there are experimental difficulties involved 
in getting an exact measurement of this parameter, the discrepancy between the theoretical and 
experimental estimates may well reflect a lower affinity of integrins when membrane bound, 
possibly due to steric constraints which will have a profound impact on this allosteric protein. 

We will therefore follow Moy and co-workers (1999) and use rj = 10 22 ^|^ for converting in- 

tegrin dissociation constants. For 77 = 10 22 ^|p we find k on ~ 2 x lCT^m 2 (J) s)^ 1 . In order to 
determine the frequency of integrin-ligand bond formation the density of open integrins needs 
to be taken into consideration. The average total (closed and open) integrin density on the cell 
surface has been estimated as p ~ 1 — 3 x 10 2 fim~ 2 [Wiseman et al. , 2004] . At least in some 
cells (e.g. platelets) more than 95% of all in tegrins are in the closed inactive conformation in 
the absence of ligand UTadokoro et al, 2003| |; the density of integrin in the open conformation 



is therefore small (p Q < 10/xm~ 2 ). Given that k + = k on p this implies, using the experimental 

estimate for rj, that k + ~ 1CT 1 — 10~ 2 s _1 such that ^ « 1, which is insufficient to drive 

integrin activation. A larger k + can be achieved if a larger fraction of integrins is in the open 
conformation already in the absence of ligand. However, this leads to ligand-independent in- 
tegrin activation. We can conclude that for experimentally determined parameters the model 
predicts that ligand-dependent integrin activation will not occur without further supporting in- 
teractions. 



A stabilizing factor is not sufficient for ligand-dependent integrin activation 

Proteins that bind and stabilize the active integrin conformation have been suggested to be 
important for integrin activation (Fig. |2j3). The Master equation derived in the previous section 
can be extended to include such factors that stabilize integrin-ligand binding. While talin is 
an excellent candidate for such a protein (and we will therefore call this factor talin in the 
following), the analysis is kept general enough such that it could be extended to any cytoplasmic 
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protein that prolongs the open active conformation by binding to the active integrin. 

Much as in the case of integrins, talins can be considered to exist in two forms of different 
activity. Cytoplasmic talin is inactivated by self-interactions and only the 'open' conformation 
can bind to the membrane (PIP2), integrins and other proteins. Therefore we can again simplify 
the model by only considering the open active form, whose concentration can again be taken to 
be constant, since any open talin bound to integrins can be expected to be rapidly replenished 
from the pool of closed talins. Given the low abundance of both open talin and integrin, com- 
plexes of the two are taken to be absent in the absence of ligand. The gains and losses of the 
state with i integrins and t talins bound can again be translated into the linear Master equation 



p it t = k_(i + 1 - t)p i+ltt + k*_ (t + l)p i+M+ i (3) 
+k + (L - i + i)pi-i, t 
+Z_(t + l)Pi,t+i + - t + l)pi, t -i 
-(k-(i -t) + k*t + k + (L + Z+(z -t) + l-t)p ijt 

with p i>t = if t > i; t, i < 0; t, i > L. Here p i)t is the probability that i (open) integrins 
and t (open) talins are bound to the ligand with L sites. p i)t = if t > i reflects the fact that 
talins can only attach to the ligand indirectly by binding to integrins. k + , Z+, Z_ refer to the 
on- and off-rates of the integrins and talins respectively. The on-rates include the density of 
integrin and talin which are taken to be constant (see above). k*_ refers to the integrin off-rate 
when talin is bound (A;* <kJ). 

For the steady state (p itt = 0) we obtain from © 





Pi,t = " 7 , x / / . x ( 4 ) 

Ei=L,t=L 
i=0,t=0 

with a = fc ~ ( ' fc ~^~| + ' +fc ~ , b = ^j— for i > t and p ijt = for % < t. The expectation 
value for the number of integrins and talins bound to a given ligand with L binding sites is 
given by < I b >= ELo^LoPht and < T b >= J2t=o t T l LoPi,t respectively. For strong 
ligand binding in the absence of talin binding (jo ij0 3> Poo) we require b > a and we recover 
the condition k + 3> k_ found in the previous section. If talin binds the condition becomes 

J2t=o PL,t ^> Po,o and thus a < b + 1. This condition can be met if either a < b or a < 1. 
Assuming that the integrin-ligand off-rate in the presence of talin is very small (fc* <C 1), a < b 
still requires k_ < k + (high integrin-ligand affinity) which we have found above to disagree 
with experimental estimates. The condition a < 1 can be met if either k- <ti k + or if k^ > k + 
and Z_ <C l+, that is integrin activation would be possible despite a low integrin-ligand affinity 
if the integrin-talin affinity were sufficient high (Fig. HJ). However, such high integrin-talin 
affinity would also lead to integrin activation in the absence of ligand binding. We can therefore 
conclude that a stabilizing factor alone is not sufficient for ligand-dependent integrin activation. 



A positive-feedback is required for ligand-dependent integrin activation 

The remaining two mechanisms in Fig. |^]both involve a positive feedback. Here, Figure 
considers a feed-back mechanism that only involves the integrin-ligand pair itself while accord- 
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ing to Figure |3) a larger network would be necessary. Both mechanisms are analysed in the 
following. 



A positive feedback that is based on the integrin-ligand pair 

A positive feedback only involving the integrin- ligand pair (Fig . I2TJ) could either be enabled by 
integrin-integrin interactions in the open form [ |Li et ah, 2003| or by a ligand-induced confor- 
mational change that leads to a higher affinity of binding. In both cases this positive feedback 
ought to be triggered once a certain (small) number of integrins is engaged in close proximity. 

Such interaction can thus be captured in the model by replacing with kL in © and setting 

f k° 

- = 71 — fc^» r » where k°_ refers to the integrin-ligand off -rate in the absence of a positive 

feed-back. k c determines the strength with which the feedback reduces the integrin-ligand off- 
rate, and K and n limit the effect of integrin-integrin interactions to local interactions. The 
model does not contain any spatial information and it is thus assumed that integrins either bind 
to pre-clustered ligand, or preferentially in the vicinity of other bound integrins. Note that the 
assumption that ligand-integrin interactions are short-lived in the absence of integrin-integrin 
interaction will lead to such preferential binding in the vicinity of bound integrins. For the 
steady state {pi = 0) we then obtain from CQ) 

(it) n}=o ( x + JJ+kt) 

F% ~ 7 l \ i Z ' ^ 

S/=0 ( l J (f ± ) rij=0 (1 + (j k +K)n) 

Given that K and n need to be chosen such that only local integrin-integrin interactions 
reduce the ligand-integrin off-rate, for sufficiently large k c © can be approximated by 




© and © are similar and only differ in the factor k c . The same argument thus applies 
such that k c 3> ^ ~ 10 2 — 10 3 is required to enable ligand-dependent integrin activation (Fig 

EJl. In the analysis of © k + was determined for the high affinity integrin-ligand interaction 
and a further ligand-induced affinity increase by 100-1000 fold is impossible. Integrin-integrin 
interactions strong enough to reduce the high affinity dissociation constant by a further factor 
of 10 2 — 10 3 are also unlikely. Thus a feedback loop only involving integrins is unlikely to drive 
ligand-dependent integrin activation. 

A positive feedback loop based on a regulatory network can enable ligand-dependent in- 
tegrin activation 

The last mechanism (Fig. Et*) to be analysed is one that considers a regulatory network that 
mediates the positive feedback. This mechanism involves a stabilizing factor whose activity 
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is increased in response to ligand binding. Such positive feedback has indeed been reported 
in the form that ligand-bound integrins as well as talins trigger an increase in PIP2, which in 
turn increases t he recruitment of talin to the membrane, and the reby the talin-integrin complex 
formation (Z + ) [ Garcia- Alvarez et ah, 2003]|Martel et ah, 2001| . The talin density on the mem- 
brane will thus be very small when integrins are inactive and increase when ligands increase the 
integrin activity. This positive feedback can be incorporated into the model as 

/ -I (t + ir (7) 

such that the talin on-rate l + depends on the number of talins in the integrin-ligand complex 
with K being the Hill constant, n the Hill coefficient and l c some proportionality factor. While 
many other formulations of such feedback are possible, this saturation form captures the likely 
talin dependence of the membrane PIP2 concentration, which will be low below a certain talin 
threshold and eventually become saturated. 

The steady state for p itt is now more difficult to derive but numerical studies show that 
such positive feedback indeed enables ligand-induced integrin activation as long as Kj L < 1 
(Fig. 13). Note that, as before in the case of integrin-integrin interactions, ligand-dependent 
integrin activation requires the proportionality factor l c to be of order 10 2 — 10 3 . In case of talin 
recruitment this is reasonable since local PIP2 production may lead to a 100-1000 fold increase 
in the local membrane talin density. A ligand-independent increase in l + , due to a signaling 
dependent increase in a factor that stabilizes the open integrin conformation, will also lead to 
integrin activation. This is likely to provide the mechanistic basis for inside-out signaling. 

We thus conclude that ligand-dependent integrin activation requires 3> k + , /_ ^> l + in 
the resting and <C k + , /_ <C l + in the active state. The first condition is met by the binding 
of the stabilizing factor which reduces and the second condition is enabled by a positive 
feedback that leads to an increases in l + in the presence of ligand. 



Conclusion 

The analysis of possible mechanisms for integrin activation suggests that a positive feedback is 
required for ligand-dependent integrin activation. Thus the ligand-integrin affinity appears to 
be too low to stabilize the active integrin conformation in the absence of a further stabilizing 
factor such as talin. To ensure that ligand-independent integrin activation by such a stabilizing 
factor alone is impossible a positive feedback that upregulates the stabilizing factor upon ligand 
binding is required. 

This conclusion strongly depends on the order of magnitude of the conversion factor 77 be- 
tween the 2- and 3-dimensional dissociation constants and further careful measurements of this 
parameter will be important. While ligand-dependent integrin activation would be possible in 
the absence of further stabilizing interactions if the theoretical rather than the experimental 
estimates for 77 were correct (that is if 77 were 1000 times smaller), available experimental infor- 
mation is very much in favour of the experimental estimate that has been used in this study. In 
addition to the argument already given (see above) further experimental observations corrobo- 
rate the notion that the physiological integrin-ligand affinity is tuned such that integrin activa- 
tion depends on supporting factors. Thus mutations in the /3-tail that disrupt the talin-integrin 
interaction markedly reduce the fraction of active integrins ITadokoro et ah, 2003| . This reduc- 
tion can be overcome by binding of an antibody that stabilizes the high affinity conformation 
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[ Tadokoro et al. , 2003| which is in agreement with the model prediction that a higher ligand- 
integrin affinity can enable talin-independent integrin activation. With integrin activation bal- 
anced on a knife-edge small changes in the affinity, density and conformational bias of integrins 
can have a large imp act on the cell's adhesiveness and motility - as is observed in experiments 
MPalecekg? al, 1997| . 

An advantage of a positive feedback regulation is the regulatory potential. Signaling pro- 
cesses can affect the positive feedback, providing a basis for inside-out signaling. It is likely 
that the observed integrin clustering is a result of such a positive feedback mechanism; a local 
increase in the density of a stabilizing factor will facilitate the formation of further integrin- 
ligand bonds in the vicinity of existing ones. Clustering would then result from a process that 
enables long-lived integrin-ligand interactions but is not of itself essential. The concept that 
ligand binding and clustering are independent processes is in agreement with the observation in 
Drosophil a that integrins can still bind to the extracellular matrix in the absence of talin but fail 
to cluster HBrown etal, 2002| . 

The feedback loop that is considered in the current model only comprises a stabilizing factor 
(such as talin) and a way of increasing the concentration of the stabilizing factor, such as by 
PIP2 formation. Many more players are known to be involved and understanding their relative 
contributions will be important. Here, especially the role of PIPKI7 will be of interest. While 
PIPKI7 produces PIP2 and thus incr eases talin recruitmen t to the membrane, PIPKI7 itself is 
recruited to the me mbrane by talin [|Di Paolo et al. , 2002| and competes with the integrin (3 
tail for talin binding [de Pereda et al. , 2005| . Integrin activation leads to FAK activation and a 
FAK-enhanced Src-mediated phosphorylation of PIPKI7 further increases its affinity for talin 
[ |Ling et al., 2002[ |Arias-Salgado et al., 2003| |. PIPKI7 therefore appears to play an important 
role both in integrin activation and the turnover of focal contacts. A combination of modeling 
and experiment, of the sort described here, is expected to shed more light on this complex 
regulatory network. 
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Figure 1 : A simple model of integrin activation. Integrins exist in either a closed (LHS) or open 
conformation (RHS); the open, high affinity form, can be stabilized by binding of extracellular ligand or 
intracellular proteins. The open conformation triggers downstream signaling and is termed active. 
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Figure 2: Four possible mechanisms for integrin activation. (A) Ligand-binding is sufficient for 
integrin activation; other factors are only important for down-stream signaling events. (B) A stabilizing 
factor (such as talin) is necessary in addition to ligand binding. (C) A ligand induced positive feedback 
(based for instance on integrin self-interaction) induces integrin activation. (D) A positive feedback 
involving further signaling molecules is necessary for integrin activation. The symbols represent simple 
variants of those in Figure ^ L refers to ligand, S to a stabilizing factor. The rate constants correspond 
to those used in the models. 
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Figure 3: Ligand engagement is not sufficient for integrin activation The integrin saturation of the 
ligand, < ^ > , is plotted against the integrin-ligand association constant K a = r 1 * 1 ; L=10. As discussed 
in the text, experimental data suggest that such that the fraction of activated integrins would be 

very small. 




Figure 4: A stabilizing factor is not sufficient for ligand-dependent integrin activation The fraction 
of ligand bound by integrin, <I £ > , is plotted against the integrin-ligand association constant K a = ^ t 
and the integrin-talin association constant K a = ^ for k*_ = 0; L=10. For experimental estimates of 
these association constants the activated fraction is small. 
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Figure 5: A positive feedback enables ligand-dependent integrin activation The fraction of ligand 
bound by integrin, <I £ > , is plotted in the absence (- -) or presence of a talin feedback regulation against 
the Hill constant K (Eq. 0> normalized by the number of binding sites L for different Hill coefficients 
[n = 1 (-); n = 5 (• • •)]. Parameters were set to L = 10, k + = 0.1, l c = 10 2 , /c_ = 5, fc* = 10" 3 , 
/_ = 1. It can be seen that as long as positive feedback is initiated upon binding of few integrins (small 
-jO, a large active integrin fraction is obtained. 
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